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Abstract. We analyze the time evolution describing a quantum source for nonin- 
teracting particles, either bosons or fermions. The growth behaviour of the particle 
number (trace of the density matrix) is investigated, leading to spectral criteria for 
sublinear or linear growth in the fermionic case, but also establishing the possibil- 
ity of exponential growth for bosons. We further study the local convergence of the 
density matrix in the long time limit and prove the semiclassical limit. 



1. Introduction 

Particle sources are an indispensable part of any scattering experiment. Nevertheless in 
the theoretical description they are mostly disregarded on the basis that a "suitable" 
wave function has been prepared. Of course, a fully realistic modelling of a particle source 
will be difficult and possibly of marginal interest. However on intermediate grounds, 
having a simple model source could be of use. The purpose of our paper is to study a, 
in a certain sense, minimal model. Surprisingly enough, at least to us, we will find that 
for bosons there is a dynamical phase transition. 

On the classical level a particle source is easily modeled and, in variation, used 
widely without further questioning. To explain the principle let us discuss particles in 
one dimension with position Xj G R and velocity Vj 6 R, j — 1,2,.... The source 
is located at the origin and turned on at time t — 0. Particles are created at times 
< ti < t% < ... . Once the j-th particle is created at time tj, it moves freely as 
Xj(t) = Vj(t — tj), t > tj. To have a Markov process it is assumed that — tj) are 
independent and exponentially distributed with rate A. Also at the moment of creation 
the velocity is distributed according to h(v )dv independently of all the other particles. 
The average density / on the one-particle phase space is then governed by 



Of course, one could imagine some statistical dependence. But then the simplicity of 
equation (11.11) is lost. In applications, the left hand side of the transport equation may 
contain further items as an external potential, a nonlinear collision operator, and the 
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like, to which the source term on the right hand side is simply added. The underlying 
reasoning for the source term still follows from the statistical assumptions stated above. 

On the basis of (j 1.1 ft one concludes that, provided f(x, v, 0) = 0, the particle number 
increases linearly in time, 

f{x,v,t)dxdv = tX I h(v)dv. (1.2) 



One can also show, that /(•,£) converges to a steady state as t — > oo. 

In the quantum case the Markov process on the many particle level is replaced 
by a quantum dynamical semigroup on Fock space. As in (II. 11) . the crucial constraint 
comes from the condition to have a closed equation on the one-particle level, a condition 
which essentially determines the model uniquely. Particles are created in the pure state 
(f) EH = L 2 (R d ) . Once created they move in Mr according to the Schrodinger equation 

ij^ = Hip. (1.3) 

The one-particle Hamiltonian if is a self-adjoint operator with domain T>(H) C TC. For 
later purpose, P pp , resp. P ac , is the spectral projection of H onto the pure point, resp. 
absolutely continuous, part of the spectrum. Primarily, we think of the free Schrodinger 
evolution, in which case H = —A (we use units such that the Planck constant h equals 
1 and the particle mass m equals 1/2). But some of our results also hold abstractly, 
in particular for the Schrodinger operator with a potential, H = —A + V(x). We use 
the shorthand H = —A. Our techniques generalize without great efforts to the case of 
particles being created in a mixed state. Hence we decided to stay with a pure state as 
minimal model. 

To be more precise, we introduce the Fock space 

oo 

f± = (BS ± H® n , (1.4) 

ra=0 

where S±TL® n denotes the (anti) symmetrized n-fold tensor product of TC with itself, i.e. 
the either bosonic (+) or fermionic (— ) n-particle subspace. The bosonic or fermionic 
creation and annihilation operators are denoted by a± (</>)* and a±(4>). They satisfy the 
(anti)commutation relations 

[a±(f),a* ± (g)} T = (f,g)lf ± , 

[a±(f),a ± (g)} T = 0=[al(f),a* ± (g)] T [ ' } 

with (•, •) the scalar product of H and [A, £]_ = [A, B] = AB — BA, [A, B] + = {A, B} = 
AB + BA for operators A, B on T±. We postulate an evolution equation of Lindblad 
type 

jA{t) = £A(t) (1.6) 

with bounded A(t) on T±, where £ is the generator of a completely positive dynamical 
semigroup. To have a closed equation on the one-particle space, £ has to be quadratic 
in a±(<f>)*, a±(4>), in other words £ has to generate a quasifree dynamical semigroup. 
Clearly £ is the sum of the hamiltonian part £ and the source part £ s . £ is obviously 
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quasifree and completely positive. If particles are created in the state G H, \\<p\\ = 1 
with rate |A| > 0, then the only possible source term has to be of the form 

C S A = |A| (2a ± (<f ) )Aa ± ( ( f ) y - a±(0)a±(0)M - Aa±(0)a±(0)*) (1.7) 

Let ujq be the initial state as density matrix on T±. Then Ut(A) = u>o(e ct A) defines the 
state at time t. Its one-particle density matrix is given by 

(g,p(t)f)=u H (a*(f)a(g)). (1.8) 

Here p(t)* = p(t), p(t) > 0, tr (p(t)) < oo and p(t) < 1 in the case of fermions. If u is 
quasifree, then ui t is also quasifree. In particular, u> t is uniquely determined by pit). 

From fll.7p . (11.81) one readily obtains the evolution equation for the one-particle 
density matrix pit), 

j t p(t) = -i[H,p(t)} + 2|A|P^ + A(i>(*) + p(t)Pt), (1.9) 

where P^ = \<f))(<f)\ is the orthogonal projection onto 0. We regard A as real parameter, 
A G M. Then A > in (11.91) is the evolution equation for bosons, while A < refers 
to fermions, and we maintain this convention throughout. (I1.9P holds for t > and 
is supplemented by the initial p(0) = po- It is the quantum analogue of the classical 
equation (II. ip . 

The goal of our paper is a detailed analyis of Equation (11. 9p . 

While precise conditions will be given in the main part, let us explain already now 
the rough overall picture which emerges from our study. We take H = H = —A, a 
"reasonable" wave function 0, and start with p(0) = 0, which corresponds to the Fock 
vaccuum. Hence ut is quasifree and gauge invariant. Basically there is a competition 
between the speed of transport through —A and the rate |A| at which new particles 
are supplied. For small |A| the Laplacian dominates. At time t the front particles have 
travelled a distance of order t from the origin. In essence, the state u> t is an incoherent 
mixture of single particle wave functions, similar to the classical set-up described by 
(1.1). If A — > — oo, i.e. per unit time a large number of fermions are created, then the 
replenishment becomes constrained because of the exclusion principle and, while the 
number of particles still increases linearly in time, the particle current should level off. 
In fact, we will show that the current vanishes in the limit of large production rate. 

On the other side for bosons, beyond some critical value A c , the operator iH Q + XP^ 
attains an isolated, nondegenerate eigenvalue satisfying (iH Q + XP^cpx = o(A)0a, 0a G 
7i, with Rea(A) > 0. Therefore in the long time limit 

p(t) = P 0A exp[2Rea(A)t]. (1.10) 

The particle number increases exponentially and a pure Bose condensate with conden- 
sate wave function A is generated. (f)\ — > (ft as A — » oo. We will prove this scenario for 
sufficiently large A. The critical regime, which in principle could be more complicated, 
remains to be explored. 

Quasifree dynamical semigroups were introduced in [51 Sj. A very readable account 
is the review by Alicki in [2] . He writes down Equations (ll.7p and (11.91) for the case of a 
general source and sink. He also discusses coupled quantum fields, for which such kind 
of equations arise in a weak coupling limit. The case of a sink only is considered in [TJ. A 
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discrete time model is studied, which however in a continuous time limit converges to a 
quasifree dynamical semigroup of the type (II. 7p with creation and annihilation operators 
interchanged. We also refer to [3J, where more recent mathematical contributions are 
listed. 

To give a brief outline: In the following two sections we properly define the solution 
to (11.91) and list the main results. In Sections H] to [6] we study the particle number, 
N(t) = tr(p(t)), when starting with an empty space po = 0. In particular we establish 
both asymptotically linear and exponential growth depending on the parameters. In 
Section [7] the convergence of the local density matrix as t — > oo is investigated. Finally 
we prove that in the semiclassical limit, |A| = 0(e) and time, space 0(e~ l ) the solution 
of equation (II. 9p converges to the solution of the classical source equation (II. ip with 
h{v)~\^v)\\ 



2. Existence of solutions 

The formal solution of equation (11.91) for t > reads 

P (t) = Sgn(A) ^(~iH+XP,)t e (iH+XP,)t _ ^ + e (-iH+xp,)t poe (iH+\p,)t^ ( 2 _-Q 

where 1 = 1^ is the identity map on 7i. Since H is selfadjoint, iH generates a strongly 
continuous unitary group on TL and, considering ±AP^, as a bounded perturbation of 
this generator, we can apply Theorem 2.1 of [6], Chapter IX., §2 to conclude that 
T = iH ± XP^ : T>(H) —>■ 7i still is the generator of a strongly continuous group which 
is norm-bounded as ||e* T || < e' A 'l. Therefore, the (semi)groups occuring in (12. ip are 
well-defined. Furthermore, the well-known formula 

e (*H + W = s _ lim ( ! _ (±iH + XP,)t \ n 



implies 



for all t E R. 

Next we have to ensure that p(t) is actually a density matrix, i.e. a positive trace 
class operator for all £ > 0. From (12. ip one obtains 

p(t) = 2|A| /* ^H+XP^sp^iH+XP^s^ + e (-iH+XP^t poe (iH+XP,)t ^ 

Jo 

in which the integral is a trace-class valued Riemann integral (recall that the semigroups 
are strongly continuous). Since po and are positive trace class operators, (12.31) implies 
the same property for p(t). In case A < 0, we have to check the fermionic property 
p(t) < 1. But if the system starts in a fermionic state, po < 1, (12.11) can be reordered to 

1 _ p{t) = e (-w+*P*)t (! _ po ) e {iH+xp+yt^ (2.5) 



which stays positive for all t. Thus the fermionic property of po is preserved in time. 
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3. Main results 

3.1. Asymptotics of the particle number 

Most basically one would like to know the number of particles, N(t), produced by the 
source. Naively one would expect N(t) to grow linearly. However, the statistics of the 
particles induces an effective attraction, respectivly repulsion, which might change such 
simplistic picture. From (12. ip one easily calculates 

j tr(p(t)) = 2|A| \\e^ H+xp ^<P\\ 2 + 2A (<f>, e^ H+xp ^ p^ H+xp ^ <p) . (3.1) 
Together with the observation that for all tp e Ti,, 

£\\ e (- iH + xp ^\\ 2 = 2A \(<p,e { - iH+xp ^)\ 2 , (3.2) 

equation (13.11) immediately implies that for any Hamiltonian H, the number of particles 
grows at least linearly in time for A > 0, and at most linearly for A < 0. The number 
of fermions is monotonically increasing, which is an easy consequence of ( 13.11) together 
with the fermionic property p < 1 an d the fact that \\e < y-' iH+xp <t>) t ( j)\\ 2 = ||e( +lH+AP "*)*0|| 2 
(their derivatives with respect to t are equal). 

Concerning the asymptotic growth of N(t), the initial density po does not change 
the qualitative behaviour, and is therefore set to p = 0. In addition to the upper bound 
N(t) < 2|A|t for the number of fermions, there is an easy characterization of those source 
states for which the number of fermions stays bounded. 

Theorem 3.1. Let A < 0. N(t) stays bounded as t —>■ oo if and only if the source state <fr 
is a finite linear combination of eigenvectors of H . In this case 

lim N(t) = N. (3.3) 

t^oo 

with N the number of different eigenvalues corresponding to these eigenvectors. 

Concerning the distinction of source states generating linear or sublinear growth of N(t) 
in the fermionic case we have the following result. 

Theorem 3.2. Let A < 0. If P pp <p = <f>, then the limit rate of particle production vanishes 
and N(t) grows sublinearly. On the other hand, if ^ 0, then N(t) increases linearly. 

In particular, for H = H , and arbitrary (j), the number of fermions increases linearly 
in time. 

For bosons, linear growth is not an upper, but a lower bound for N(t). But also 
in this case one can characterize explicitely a class of source states which yield linear 
growth for small A. 

Theorem 3.3. Let (p E Ti,, \\(f>\\ = 1 and assume that r := J °° | ^(p,e~ tHt (p^ \dt < oo. Then 
for all A < or < A < r -1 , -^N(t) converges to a non-zero limit as t — ► oo, which 
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satisfies 



2|A| < lim ±N(t) < 2|A ' 2 (A > 0), 
t^oo dt (1 — At) 



•>l\l ,1 (3 ' 4) 

< lim -j-N(t) < 2 1 A | (A<0). 



(1 - At) 2 t->oo dt 

The quantity (0, e~ lHt cf)) is the overlap between the source state at time and at time 
t. Thus r should be regarded as a measure for how long it takes the time evolution to 
transport an emitted particle away from its source. In this context, |A| < r _1 means, 
that the strength of the source is smaller than the "transport capacity" of the time 
evolution. Therefore, the emitted particles hardly influence each other, and in essence 
the bosonic or fermionic character does not show in the evolution of the density matrix. 
The limit particle production is constant as in the classical case. An example for such 
in the physically relevant case H = L 2 (IR 3 ), H = —A is G L 1 fl L 2 (M 3 ), since 
lHt (j)^ has a t -3 / 2 decay, as can be deduced from the free propagator 

(e- iHot <f)) (x) = (4mt)~% I e^^4>{y)&y (3.5) 

for G L 1 fl L 2 (]R d ) (compare [llj, Chapter IX. 7). For fermions, this is also an example 
where one has an explicit estimate for the speed of particle production. However, the 
form of the lower bound in (13.41) already suggests that the limit growth rate might 
decrease to as A — > — oo. In fact, we have 

Theorem 3.4. For G V(H), we have 

lim lim f ^N x {t) J = 0. (3.6) 



A^-oot^oo \dt 

For bosons, the behavior at large values of A is qualitatively very different. In this 
regime, exponential growth occurs for arbitrary choices of the Hamiltonian H and the 
source state 0. 

Theorem 3.5. For sufficiently large A > the operator iH + XP^ has an eigenvalue a(X) 
with positive real part, and liniA^oo A _1 a(A) = 1. 

If such an eigenvalue a (A) with normalized eigenvector if) exists, the number of 
bosons can be estimated by 

tr(p(0)>(^pW)>e 2ReQ(A)t -l. (3.7) 

Hence Theorem 13.51 implies the existence of a critical strength A c > such that the 
number of particles grows exponentially in time for all A > A c . For source states <fi as in 
Theorem 13.31 our result implies a dynamical phase transition from linear (0 < A < r _1 ) 
to exponential (A > A c ) growth. The general picture is more complicated, however. If H 
has point spectrum, it is obvious that choosing as an eigenvector of H will generate 
exponential growth for all A > 0: N(t) = e 2Xt — 1. But also for H = H Q there are source 
states which generate exponential growth for all A > 0. 
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Theorem 3.6. Let e V(H ) = H 2 (lR d ) with Fourier transform 

kp) = |^|-V. H (^^Z3Z:. (3.8) 

T/ien iP + ^P/> has an eigenvalue a (A) witt positive real part for all A > 0. 

In the proof of this theorem, the eigenvalue will not be computed explicitly. With 
a view towards the semiclassical limit, one can infer that Rea(A) = o(A) has to hold. 
In this limit one considers a source with activity A = ce on a time scale e~ l t and thus 
the exponent in (13.71) . 

2Rea(ce)e~ 1 t, (3.9) 

has to vanish in the limit e — > so to yield the linear classical growth behaviour. 

In contrast to Theorem 13.31 the source state (13.81) has an overlap decay as 1/Vi, 
thus r = oo. As a consequence, the emitted particles stay close to the source for a long 
time and, being bosonic, they pull further particles out of the source, which leads to an 
exponential growth of the particle number regardless of how small A. 

3.2. Convergence of p{t) 

Having studied N(t) one may ask whether p(t) has a limit as t —>■ oo. Since mostly 
N(t) — > oo for large t, the natural notion is to study the local limit of p(t). Let Q C M. d 
be a bounded region and let Pq denote the orthogonal projection of Tt onto the subspace 
L 2 (Q). We consider the number of particles in Q, i.e. 

N n (t)=ti(P n p(t)P n ). (3.10) 

If this quantity stays bounded as t — > oo, one can use that Pnp{t)Pn has a positive time 
derivative (at least for p = 0) to infer that the restricted density matrix P^p{t)PQ has 
even a trace class limit. 

A first result shows that for fermions and H = Hq this limit, as far as it exists, does 
not depend on the initial condition p . 

Theorem 3.7. For Q C M d , with finite Lebesgue measure |0| < oo, A < and po an 
arbitrary trace class operator, it holds 

^p ne (-iHo+XP,)t poe (iH 0+ XP,)t pn ^^ Q asf _, 00- ( 311 ) 

For all remaining results, we return to po = again. The next result is rather 
obvious: If particles are generated at most linearly in time, and if they move away from 
a certain region Q in space within finite time, the number of particles in Q will not 
diverge: 

Theorem 3.8. Let us choose and A such that N(t) has a linear bound and assume that 

POO 

/ \\P n e- iHt <j)\\dt < oo. (3.12) 
Then tr (P^p(t)Pn) approaches a finite limit as t — ► oo. 
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In particular, for H = H this theorem applies to f2 with \Q\ < oo for all 4> G 
L 1 fl L 2 (R 3 ) and all A < r _1 : (Sub)linearity follows from Theorem 13.3} and (I3.12p is 
deduced from the explicit form of the integral kernel (13.51) . 

The next result is characteristic for the behaviour of a fermionic particle source. 
Consider the fact that in a region Q of finite measure, there are only finitely many 
states with kinetic energy below a certain bound (cf. [9], p. 27). If the source state has 
no high-energy contributions, it should only be able to charge a finite number of states, 
so that the particle number in fl stays finite. 

Theorem 3.9. Let H = H , <p G TC, with Fourier transform <f> supported in K C M. d 
where \K\ < oo, and let fl C M. d with \Q\ < oo. Then, for A < ; 

lim tr(P n p{t)P Q ) < (2n)- d \n\\K\ < oo. (3.13) 

The upper bound corresponds to the phase space volume of fermions, i. e. to one 
fermion per unit cell. 

3.3. Semiclassical limit 

To pass from the Hilbert space formulation of quantum mechanics to phase space, one 
convenient tool is the Wigner transform. One defines 

W[k](x, P ) = (27r)- d J ' b-^k (x+ |,x - |) dy, (3.14) 

where k G L 2 (M. d x M.y) is the integral kernel of a Hilbert-Schmidt operator. In par- 
ticular, since all trace-class operators are Hilbert-Schmidt, we can consider the Wigner 
transform for any density matrix. Denoting, as usual, the semiclassical parameter by e, 
the semiclassical limit corresponds to the small e behaviour of 

f'(X, P, T) = e~ d W [p e (e _1 T)] (e' x X, P) , (X, P, T)el rf xR d x R+ (3.15) 

where p e {t) is the solution of 

j t p e (t) = -i [H , p e (t)] + 2|c|eP^ + ce (P^it) + p e (t)P») (3.16) 

for all t > 0, with the convention Hq = — A/2. The initial densities p e Q are chosen such 
that ||polltr is bounded uniformly in e > and 

e- d W[pl](e- 1 X,P)^g(X,P) (3.17) 

for some distribution g G V (M| x R d p ). The equations fl3~T5|) - (l3~T?]) describe a quan- 
tum particle source producing particles on a microscopic scale (x,p,t) correlated to the 
macroscopic scale (X, P, T) by (X, P, T) = (ex,p, et), the standard semiclassical scaling, 
see e.g. [S]. To have a bounded rate |c| on the macroscopic time scale, we have set A = ce. 
In this case, we have the following convergence result for the phase space density f e . 

Theorem 3.10. For all g G V (R d x R d ) , G H, \\(j)\\ = 1, and for all T > 0, the limit 
hmr(X,P,T) = f°(X,P,T)=g(X-PT,P) + 2\c\[ 5(X - P S )|0(P)| 2 ds (3.18) 
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holds in the topology of V (M.^ x M. p ) . This limit solves 

(X,P,T) + P ■ V x f(X,P,T) = 2\c\6(X)\4>(P)\\ 

f(X,P,0)=g(X,P) 
in the sense of distributions on phase space. 

(13.191) is the weak form of equation (II. ip with a source term defined through the 
semiclassical limit of <fi. 

A natural step would be to include an external potential varying on the macroscopic 
scale. We leave this as an open problem. 



4. Particle production in the fermionic case 

4.1. Proof of Theorem O 

Proof. Since we are in the case po = 0, p(t) is different iable in trace class with a positive 
operator as derivative, 

jp(t) = 2\X\e^ H+xp ^P^ lH+xp ^. (4.1) 

Thus \\p(t) — p(s)\\tr = I tr(p(t)) — tr(p(s))|, so if the particle number stays bounded and 
tr(p(t)) approaches a finite limit from below, p(t) converges in trace class as t tends 
to infinity. As a fermionic density matrix, the limit p^ obeys < p^ < 1 and it stays 
invariant under the time evolution given by (12. ip . 

Poq = 1 _ ^-iH+XP^t^iH+XP^t + e (-iH+XP^)t p ^ e (iH+XP^t (A^) 

Furthermore, since tr(p(£)) converges to a finite limit, its (monotonously decreasing) 
time derivative tends to zero, and therefore, 

(0, (1 - Poo )<P) = hm (<p, (1 - p(t))<P) = hm ||e^ +AP ^VH 2 

t— >oo t— >oo 

d ( 4 - 3 ) 

= lim He^+^Vf = lim(2|A|)- 1 ^tr(p(t)) = 0. 

t— >oo t— >oo dt 

Since 1 — p^ is positive, = (4>, (1 — Pco)<P) can only hold if (1 — poo)4> — 0. Together 
with (|4T2|) this yields: 

I. ( e im pooe -iH^ = \ X \e im P^- lH+xp ^(l - Poo )e( lH+xp ^e- iHt 

+ \X\e iHt e^ iH+xp ^(l - Poo )e^ +AP ^e- im ( 44 ) 
= \\\e lHt (P^l - Poo ) + (1 - Poo )P ) e~ iHt = 0. 

Therefore, we have for t G R 

e tHt Poo e- mt = Poo . (4.5) 

Now consider the eigenspace V of p^ corresponding to the eigenvalue 1. V has a positive, 
but finite dimension because G V and p^ is trace class, and due to (14. 5p . V is invariant 
under the action of the group e~ lHt . So since V is finite-dimensional, we have V C T>(H) 
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and V has a orthonormal basis of eigenvectors of H, and can be written as a linear 
combination of those. 

Conversely, assume that (f> = Yln=i c n' l l ) n with \\ip n \\ = 1, Hip n = u) n ip n , with all 
c n 7^ and all u) n different. —iH + XP^ is a bounded operator on W := span (-01, ...,tp n ), 
so e^~ lH+xp ^ t (j) can be calculated by the power series, to see that 



p(t) =2|A| f ' eHH+XP^sp^iH+XP^s^ 

Jo 



(4.6) 



contains only states from W, which, together with the fermionic property, yields the 
bound tr(p(t)) < N. Thus p^ exists and, by (14. 5p commutes with H on W, so that they 
have a common basis of eigenvectors, which necessarily means 



N 



Poo = J ^2b n \ip n )(ip n \. (4.7) 

n=l 

Now it only remains to check that p^ = l\y, which follows from 

N N 

^ Cntpn = (fi = Poo<P = ^ C n b n tp n . (4.8) 
n=l n=l 

and the fact that all c n ^ 0. So tr(poo) = dim(W) = N. □ 

4.2. Proof of Theorem O 

Proof. For a state <fi G 7i pp we have to show that 

lim 4- tr(p(t)) = lim 2|A| ||e ( - iH+AP ^V|| 2 = 0. (4.9) 
To this end, expand in a basis of eigenvectors of H, 



J2 C n*Pn (4.10) 



n=l 



with \\ip n \\ = 1, Hi/j n = u) n ipn- Since 

^iH+XP^t^-iH+XP^t > g 

J (4.H) 

^5. (iff+AP )t (-iif+AP )t = _2|A| e (iif+AP * )t P,e ( " !ff+AP * ) * < 

— ' 

polarization implies that e (*^+ AP </>)* e (-*^+ AP </>)* converges weakly to a bounded, positive, 
selfadjoint operator A with 

A = ^iH+xp^t ^{-iH+xp^t _ (4-12) 
As the ip n are eigenvectors of if, 

(V„, e- 4m Ae 4H Vn> = (V„, ^n) • (4.13) 
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The last two equations prove 

0=j t (TPn,e- iHt Ae iHt ip n ) 

= A e*"(P+A + AP^)e im if) n ) 
= \(if) n ,{P <j> A + AP (t> )if) n ). 

Taking the sum over all n, one obtains (0, A<p) = 0, (even A<f> = 0). By the definition of 
A, this means 

lim ||e ( -^ +AP ^Vf = Jim ((p, e (iH+XP^t e (-iH+XP^t^ = ^ A ^ = 0> ^ ^ 

which proves the first assertion. 

For the proof of the second part of the theorem, assume that there is a source state 
with P a c0 7^ and a A < such that the number of particles grows sublinearly in 
time, i.e. 

h(t) = ||e(-^ +AP ^V|| = ||e( iP+AP ^V|| ^ (t -> oo). (4.16) 
By the same arguments as in (14. lip , we have the weak convergence 

^iH+XP^t^-iH+XP^t _^ A (4.17) 

with A having the invariance property f)4.12p . and, by assumption, 

A<p = 0. (4.18) 

Following the proof of Theorem 13.11 one can deduce from (14.121) and A(f> = that 

A = e im Ae~ iHt (4.19) 

for all t. Thus A and H commute, and so do A and all operators given by the functional 
calculus for H, 

Af(H) = f(H)A, (4.20) 

for all bounded Borel functions / on M. We will now apply this fact to a certain choice 
of /. For if; = P ac (p ^ 0, we have p^(dE) = p(E)dE with a positive L X (]R) function p. 
Defining the Borel function 

fc = Hp < c} ■ 1 (R \ supp(/i pp )) • 1 (R \ supp(/i sing )) (4.21) 

for c > 0, we have 

iP c = f c {H)<f> ^if) (c -> oo) (4.22) 

and 

Aif) c = Af c (H)<P = f c (H)A<P = 0. (4.23) 
Thus for c fixed large enough, if> c ^ and 

lim 1 1 eC-"^^*) Veil 2 = (ipe, Aif) c ) = 0. (4.24) 

t— >oo 

Considering the semigroup as a perturbation of the unitary group, we can write 

j-w+XP+yt^ = e -mt^ c + x f t e (-^ + AP,) (i - s ) ^ e - iP > c ) ds. (4.25) 
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By the choice of ip c , the scalar product reads 

<0,e-^ c )= / l{p<c}p(E)e iEs dE, (4.26) 

and therefore is an L 2 (IR S ) function with norm not larger than \/2nc. Choosing e > 
and T > such, that 



| (<j),e~ iHs ip c ) \ 2 ds < e (4.27) 
we can split the integral in two parts to see for all t > T 

|| e (-iH+AP^ c || > ||^|| _ | A | j T h[t _ s) |(0 )e -^ c >| ds 

° , t (4-28) 
-|A| sup / |(x,e ( - iH+A ^ )( '- s V)<0,e-^c)|ds. 

Now actually both scalar products in the second integral are L 2 functions, since 

— \\e^ H+xp ^ s x\\ 2 = 2\\( ( p,e^ H+xp ^ s X )\ 2 (4-29) 
ds 

and thus 

e <^). x) f ds= \\x\\ 2 -^\\f" W '>'x\\\\\£ (4 30) 

i) 2|A| |2A| 

Therefore an application of the Cauchy-Schwarz inequality to (14.281) shows 



c 



{-iH+XP^t 



T 



ipc\\>Hc\\-\\\ / h{t-s)\(<f>,e- iHs ip c )\ds-\ 1 -^-. (4.31) 



o 



Ale 



For fixed T, the integral in the last line tends to zero as t — > oo by dominated convergence 
and the assumption that h(t) — > 0. Since one can take e arbitrarily small, this would 
imply 

lim ||e ( -^ +AP ^V c || > ||^ c || > 0, (4.32) 

contradicting (I4.24p . Thus sublinear growth is not possible for source states with 
P ac 0^O. ' ' □ 



5. Linear growth 

5.1. Proof of Theorem S31 

Proof. For t > define h[t) = ||e ( - _t ^ +AP ^^0||. Considering the nonunitary time evolu- 
tion as a perturbation of the unitary group, one can write 

e (-^+AP,)i = e -^ + A /" ^-iH+XP^t-s)^ <0 )e -^0) ds. (5.1) 

JO 

For < A < r _1 , the monotonicity of h(t) implies 

h(t) < 1 + XrhCt) 

(5 2) 

h(t) < (1 - Xt)-\ 1 ' ; 
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so lim^oo h(t) exists and is bounded by (1 — Ar) 1 . For A < 0, lim^oo h{t) = exists 
since h(t) decreases monotonically, and dominated convergence implies the estimate 

hoo>l- lAlr/ioo 

^^(l-Ar)- 1 . (5 ' 3) 

Setting p = in (13.1 p and using the unitarity of e iHt , this implies the existence of 
lim^oo 4:N(t) and the estimates 

, lim ~T f N( ( t ) - 7T^T^2 ( A >°) ( 5 - 4 ) 
at (1 - Ar) 

and 

lim -^iV(*) > -^^2 (A<0). (5.5) 
t^oo dt (1 - Ar) 



The other bounds follow from (13. ip and (13.21) . 

5.2. Proof of Theorem S31 

For A < 0, G V(H), one has 

d 



□ 



and thus 



Therefore, 



so that 



e (iH-XP,)t e XP, t(f) = le Xt R ^ (56) 

(JjL 

% f e Xs ^- iH+xp ^-^H(t)ds. (5.7) 
Jo 

^-w+xp^u < e -\x\t + _L ^ ( f _ OQ ) j (5. 8 ) 

A A 

£ ISf- (5 - 9) 

which is exactly the |A| _1 suggested by (13.41) . 

To see why this theorem does not hold for all <fi T>(H), consider the Hilbert space 
TC = L 2 (Ma,) with the multiplication operator H = x and the source state with 
4>{x) = , 1 In this case, the simple form of the overlap ((j),e~ lHt (f)) = e~'*' allows 

for an explicit representation 

/ (X-l+ix)t _ i 
(-iH+XP.)t,\ / „\ _ -ixt \ I x£_ 



(e ( - iH+AP ^V) (s) = e-"V(z) ( Al Y31 — + 1 ) ( 5 - 10 ) 



for all t > 0. As t — > oo, for A < r = 1, e strongly converges to the limit function 

• (5.H) 
A — 1 + ix 



with norm 1 / \J\ — A, yielding the limit 

lim lim ( — A^.„, , 

A^-oot^oo yctt J A^-oo |1 — A| 



lim lim I ^-A" A (t) J = lim t-^L = 2 , (5.12) 
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which is a positive saturation value for the particle flux. Also note that for the critical 
coupling constant A = 1, 

(e^+^V) (x) = e' ixt <\>{x) (l + 6 "*~ (5.13) 

with 

\\e { - iH+p ^(t)\\ 2 = l + 2t (5.14) 

and thus Ni (t) = t 2 + 1. So in this case, the transition region implied by the theorems 
13.31 and 13.51 really consists only of the critical point A = r = 1, with quadratical growth 
of the particle number at the transition point itself, and exponential growth behaviour 
of the form e^ -1 )* for larger values of A. 



6. Exponential growth 



6.1. Proof of Theorem E5] 

Proof. It suffices to show that for any choice H, <p and < e < | the operator 

™+P* (6-1) 

has an eigenvalue inside the open ball B e (l) for sufficiently large A. The idea of the proof 
is, that, given a closed operator A on a Banach space B, and a closed curve T in the 
resolvent set of A such that it separates the spectrum of A in two parts o\ inside and 
o"2 outside the curve, the operator defined by 

P = Pa.y = ^~ f(z- A^dz (6.2) 



r 



2vrz 

is a (non-orthogonal) projection with the properties 

PA C AP with a (A\PB) = 0i, 

(1 — P)A c A(l — P) with a(A\(l-P)B) = a 2 . (f> " >) 

Now an application of the Neumann series shows that if B is a perturbation of A with 
|| -B (A — ^) _1 || < 1 for all z G T, one has 

dim(P Air £) = dim(P A+B:F B). (6.4) 

In particular, if the spectrum of A enclosed by T, <T\{A), consists of one eigenvalue of 
single multiplicity, then so does cr\{A + B). The details are given in Chapter IV, §3 of 
|6j. 

Setting T = dB e (l), we apply this theorem twice, first with B = He, a suitable 
"low energy" subspace of H, and then with B = TC to remove the energy cutoff again. 

Since H is selfadjoint, one can define the spectral projection Q E = 1[_ E;E }(H) by 
the functional calculus for any E > 0, and write 4>e = Qe4>- Since Q E — > In strongly 
as E — ■> oo, there exists an E > with — (p E \\ < e/8. For this choice of -E 1 , take 



1~~l E = Q E H and A > AE/e. Writing P ( f >E = \4>e)(4>e\, we first consider the operator 

iH , N 

— + 4*e (6-5) 
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restricted to He- P<f> E has exactly one eigenvalue in B e (l), which is ||0_e|| 2 £ [1 — e/4, 1] , 
and thus we have the estimate 

max \\iH (P^ - z)- 1 /X\\ < j ■ 1 < ~. (6.6) 

Therefore, the above-mentioned perturbation result applies, and the operator from (16. 5|) 
has one single eigenvalue in B e (l). Furthermore, one can estimate its resolvent by the 
Neumann series to obtain 



max II (iH/X + P <t>E - z) 1 II < max II (1 He + iH (P^ E - z) 1 /A) 1 II 

• II (p*. - -r 1 II < 6 - 7 > 

< (1-1/3)" 1 - 4/(36) = 2/e 

on the restricted space He- On the orthogonal complement He, we have iH/ A + P ( / >E = 
iH/X, so this operator has purely imaginary spectrum, and its resolvent can easily be 
estimated by 

max || (iH/X + P fe - z)~ l \\ = max II (iH/X - z)' 1 \\ < —— < 2. (6.8) 

zgr ■ zer 1 — e 

Since the operator iH/ X + P$ E on H is decomposed by the pair He, H e as described in 
[6], Chapter III., §5.6, one can simply combine the results concerning the subspaces, to 
obtain that the operator on the whole space still has one eigenvalue in B e (l) and obeys 
the estimate 

max || (iH/X + P^ E - z)' 1 \\ < 2/e. (6.9) 

Now we are ready for the second application of the perturbation result. As ||P</> — Pc/> E || < 
2||0 — (/>e\\ < e/4, we have 

max || (P, - P+ a ) (iH/X + P fe - z)- 1 \\ < - ■ - = \. (6.10) 
zer 4 e L 

As a consequence, iH/X + P^ E + (P^ — Pfe) = iH /X + P^ has one eigenvalue of muliti- 
plicity one in B e (l). □ 

6.2. Proof of Theorem S3] 

Proof. The polynomial 

p(z) := 2 4 + 2iz 3 + (Xi - 2)z 2 - (2A + z)z - (3/2)Ai (6.11) 

has a root with positive imaginary part for all A > 0. To show this let Zk, k — 1, 4 be 
its roots, counting mulitplicities. For all 2 with p(z) ^ 0, the product rule implies 

4 



p(z)p'(z) p'(z) •sr-^ 1 



IpO)I 2 pO) Is* - * 2 * 



Et^ < 6 - 12 > 



so if all Imzfc < 0, the right side, and therefore also p(z)p'(z) would never have a positive 
imaginary part for zsl. But 



Im 



(p(^Xj2)p'(^X/2)) = A/2 + o(A), (A - 0), (6.13) 
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so p must have a root with positive imaginary part for small positive A. Furthermore, p 
does not have a real root for any positive A, because one can easily check that the real 
and imaginary part of 

pit) ■ {f - 2it 2 - 2t + i) = (t 6 + 1) ■ t + 3A/2 + \{2it b + it 3 + 2it) /2 (6.14) 

cannot equal zero for the same i 6 I. Since the roots of p depend continuously on A, 

this means that there is a z with lmz > and p(zq) = for all A > 0. 

The vector ip e T*(Hq) which we will prove to be a (non-normalized) eigenvector is 

</>:= (#o-<TV (6.15) 

Taking Cr = {\z\ = R,lmz > 0} U (—R,R) as contour of integration, one has by 
Cauchy's Integral theorem and by the fact that Im z > 



lim — / — — — — ^rdz 

K-oc 2vr J Cr {z* + 1) (z 2 - 4) 

. iz\ - 2z - 3i/2 



(6.16) 



Zq + 2izq — 2z 2 — iz 



i 

= "A' 

where we have used p(z ) = in the last line. Therefore, ip is an eigenvector: 

(2 \ — 1 2 
#o-y) <t> + \{<f>,ip)<f> = i<t> + i^'ip-i<t> = a{X)ip. (6.17) 

Since H has no eigenvectors, we have P^ijj ^ 0, and thus 

|2 Rea(A) = Re (?/>, (iH + AP ) ip) = A| (<f>,ip) | 2 > 0. (6.18) 

□ 



7. Restricted limits 

7.1. Independence (proof of Theorem 13.71) 

For the proof of Theorem I3.7[ we start with two lemmas: 

Lemma 7.1. For A < 0, an arbitrary normalized source state, 

e (-iH +XP^)t _^ q /j-Q 

in the weak operator topology as t — > oo. 

Proof. Let r/>, x W. We have to show that ^0, Q{- lH o+ xp 4>) t ^ — > o as i — ► oo. To do so, 
we write the scalar product as 

e(-^ +AP ^* X > = (V, e-^°*x) + A f e^-ty) <0, e(-^ +AP ^x) ds (7.2) 

Jo 
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The first term on the right-hand side tends to as t — > oo by the Riemann-Lebesgue 
Lemma (cf. |9J|, p. 18, Example 1.1, note that TC^Hq) = H, i.e. H has only absolutely 
continuous spectrum). The second term will only be analyzed for all if) from a dense 
subset of Ti. Since (if), e (- lH o+^P^)t^ depends continuously on if) uniformly in t > 0, this 
is enough to prove the assertion. 

According to Theorem 1.3 in [9] (p. 20), there is a subset JC(H ) which is dense in 
7~La.c.(Ho) = TC such that 

(if),e mt (f>) eL 2 (R t ) (7.3) 

for all if) G K{Hq). Thus, by (I4.30p . the integral on the right side of (17. 2p can be estimated 
by the convolution of two L 2 (M) functions, and therefore tends to as t — > oo. □ 

Lemma 7.2. For a set Q C M. d with finite Lebesgue measure < oo, and A < 
in strong operator topology as t — > oo. 

Proof. If q > i,q > 2, one has l u (x) G L q (M%) and (j+i) * e L q {R d ), and by 
Theorem XI. 20 in [T2], (p. 47), 

P n (#o + i)" 1 = (^T^ + eJ q = {Ae B(H) : tr < oo.} (7.5) 

Thus Pn(Po + *) _1 i s compact. First, take if) G T>(H ) = H 2 (M. d ) and apply the fact 
that a strongly continuous group commutes with its generator, 

p ne (-iff +AP^ = p n ( Fo + iyi(H + iXP^ + i - zAP )e ( - iHo+A ^ ) V 

= P n (H + z)-V-^°+ Ap ^(# + zAP^ + i)if) (7.6) 

- zAP n (#o + i^P^-^+^tp. 

Both terms on the right-hand side tend to in || • ||^ as t — > oo. The first one because 
of Lemma [7.11 and the fact that the compact operator Pq (Hq + i) -1 maps weakly con- 
vergent sequences to convergent ones, and the second one directly by Lemma \7. II Since 
T>(H ) is dense in TC, this already implies the claim by continuity. □ 

Now one can easily show that the initial density po does not contribute to the value 
of lim t _ )00 tr(PQp F (t)PQ) in the fermionic case. 

Proof. (Theorem \3. 7| ) Since s — lim^oo p n e^~ lHo+xp ^ t = and p G J\ (the trace class), 
Lemma 3.1 from [8 J applies, and one has: 

||p ne HHo+AP )t poe (iHo+AP )tp n || tr < HP^-^+^Volltr ~> (7.7) 

as t — > oo. □ 
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7.2. Existence (proof of Theorems [Ml \SM 



The main idea for the proofs of Theorems 13.81 and 13.91 is the formula 

j t ti(P nP (t)P n ) = 2|A|||P^- iP °+ AP ^0|| 2 , (7.8) 
which is obtained from (12. ip when setting po — 0. 

Proof. (Theorem \3.8\) The assumption of sublinear particle production (which is always 
the case for fermions) implies by (13.11) and (13. 2p . that 



{ -iH 0+ \P,) S( p\ |2 ds < ^ (7>g) 



Furthermore, 



Pae' iHot (j) + A / P n e- iHoit - s] <j) (<f), e (- iP ° +AP ^» ds 
Jo 

< ||^e-^°Vll + |A| / \\Pne- iHoit - s) <f)\\\ (0, e ( - lPo+AP ^0) |ds 
Jo 



(7.10) 



By assumption (13TT21 . \\P n e- iHot <f)\\ e L 1 n Thus, the first term in dTTTUD obvi- 

ously is in L 2 (M.f ), and by (17. 9p . the second one is a convolution of an L 1 (M + ) function 
with an L 2 (R + ) funct ion, and therefore also contained in L 2 (M.f). By (17.81) . this proves 
the assertion. □ 

Proof. (Theorem \3.9\) Let K : = supp(0) be the support of 4>. Since Fourier transform 
diagonalizes H , 



J-iHo+XPjtf = e -iH t + Xe -iH t / e iff o *0 ^ e (-iH +\Pt)s^ dg ( ? X1 ) 

Jo 

implies supp (j r e^~' lHo+xp ' l '^ t (f)) C K for all times t > 0. This means 

p^i-iHo+xp^t^ = p^i A ,(_iV)e(- iHo+AP ^V, (7.12) 
where PqI^ (— iV) is a Hilbert-Schmidt operator with corresponding norm 

||PnM-<V)|U < (2^)^ /2 (l^ll^l) 172 (7.13) 
by theorem XI.20 in [12] . p. 47. Applying (I4.30p to a singular value decomposition of 
-Pol_ft:(— iV) one obtains 

/•oo 

hm tr(P nPjP (t)P n ) = 2|A| / WPne^+^fdt 

||P f ,U(-*V)e(-^ +AP ^Virdt (7.14) 

< \\P n lK(-iV)\\% 

< (2iry d \tt\\K\. 

□ 
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12) 



(8.3) 



8. Semiclassical limit (proof of Theorem 13.101) 

To obtain a semiclassical limit, it is important to observe that the distinction between 
the bosonic and fermionic character of the semigroups disappears in the limit e — > 0. 

Lemma 8.1. For any 0, if) G TL, ||0|| = 1 one has for all t,c £ R: 

lim || exp(-iH t/e + cPM - exp(-iH t/e)ip\\ = 0. (8.1) 

e— >0 

Proof. One can write this difference as 

exp(—iH t/e + cP^fjip — exp(—iH t/e)if) 

— c ((f), exp(— iH s/e + cP^s)ip) exp(— iH (t — s)/e)<f)ds. 
Jo 

For the scalar product in the integral of (18.21) . one has 

(4>,exp(-iH s/e + cP$s)ip) = (4>,exp(-iH s/e)ip) 

+ c (</>,exp(— iH r/e + cP^ip) (0, exp(— iH (s — r)/e)(f>) dr. 
Jo 

The first term in (I8.3P tends to zero for all s ^ as e — ► by the Riemann-Lebesgue 
Lemma. By the same argument, the integrand of the second term tends to zero almost 
everywhere, and it is bounded by e' cs L By dominated convergence this implies 

lim ((f), exp(-iH s/e + cP^s)^) = 0. (8.4) 

for all s 7^ 0. But dominated convergence also applies to (18.21) and the assertion is 
proven. □ 

To show Theorem 13.101 it is useful to have the inverse for the Wigner transform, 
which is given by the Weyl quantization. For a function a e L 2 (M. d x R^J one can define 
an operator on TL by 

(Op[a]V0(z) = / / a(^-,p)e^-y^(y)dydp (8.5) 

( 27r ) Jr* J** 2 

Up to a factor, this is a unitary map from L 2 (WL d x IR d ) to the space Ji of Hilbert- 
Schmidt operators on TL (cf. equation (3.3) in |8J): 

/ a(x,p)b{x,p)dxdp = [2n) d tr (Op [a]* Op[6]) . (8.6) 

With this definition, one has for all Hilbert-Schmidt operators k, 

Op [W[k]] = j^y d K, (8.7) 

where we identified operators and kernels. This is the main tool for the proof of Theo- 
rem innni 



20 



M. Butz and H. Spohn 



Proof. Let 9 € V (R d x x R d p ) = (R d x x R d P ) be a test function. We have to establish 
the existence of 

lim / f e (X,P,T)9(X,P)dXdP. (8.8) 
By (18.61) one can transform this integral to 

/ f e (X,P,T)9(X,P)dXdP = [ e~ d W[p e (T/e)](X/e,P)9(X,P)dXdP 

= / W[p e (T/e)](x,p)9(ex,p)dxdp 

JM? d 

= tr (p e (T/e) Op[0 e ]), 

where Op [a] = Op [a]* is used, and 9 e (x,p) = 9(ex,p). 
By (12. ip . the following expression holds, 

p e (T/e) = sgn(c) (exp(-iH T/e + cP^T) exp{iH Q T/e + cP^T) - 1) 
+ exp(-iH T / e + cP^T)p e exp(iH T/e + cP^T). 



1.9) 



To keep notation simple, we first assume p Q = 0. 

The estimate || (\ip)(4>\ - \x)(x\) ||tr < ll^-xll • 
Lemma 18.11 and the dominated convergence theorem, 

r T 

||p e (T/e) -2\c\ / e-^^P^^dSWt, 
Jo 
r T 

<2\c\ 



(8.10) 

x||) yields, together with 



-iH S/ep e iH S/e 

< 2\c\ (e |cT| + 1) f || e -^ 5 / £+cP ^0 - e-^I^WdS ^ (e -> 0) 
Jo 

Furthermore, by the theorem of Calderon-Vaillancourt (Theorem 2.8.1 in [7j), the oper- 
ators Op[0 e ] are uniformly bounded in B(TC) as e tends to zero, and thus, by the linearity 
and cyclicity of the trace, 



lim 



tr (p e (T/e) Op[9 e ]) - 2\c\ [ tr ( e iH ° s/e V [9 e }e- iHoSlt P^) dS 

Jo 



0. 



(8.12) 



The Heisenberg time evolution applied to Op[# e ] can be carried over to phase space by 
defining 

r} £ s (x, p):=6 e (x + pS/e, p) = 9(ex + pS, p) (8.13) 

so that (cf. [8]) 

e iH S/e Op pe^-iH S/e = Qp^]. (8.14) 

Now Lemma 3.2 from [5] applies to rjg, stating that the operators Op[?y|] e>0 are uni- 
formly bounded in B(TC) and converge strongly to Op[?7g] as e tends to zero. Thus also 
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tr (Op [f?s]-P</>) converges and, by dominated convergence, the desired limit exists, 
Urn/" f e (X,P,T)6(X,P)dXdP = lim2|c| / tr (Op[^]P ) d5 

~oy R2d j Q (gi5) 

= 2|c| / tr (0 V {t? s ]P$) dS. 
Jo 

By (18.131) . r)g does not depend on X. Op[rjg] is a multiplication operator in momentum 
space and thus 

2\c\ [ tr (Op^Pj dS = 2|c| / / 0(P£, P)|0(P)| 2 dS. (8.16) 

At this point of the proof, one should remark that rj s is no longer in 
so Op[r?°] is not given by our definition of the Weyl quantization (18.51) . But it can be 
obtained as the quantization of a C°° (W 1 x M d ) symbol with all derivatives bounded. 
This construction is described, for example, in Chapter 2 of [7|. 

Next, we allow a general initial condition, i.e. p e is an arbitrary trace-norm bounded 
sequence of density matrices such that (13.171) holds for some distribution g on macro- 
scopic phase space. Concentrating on the last term in (18.101) . one has to evaluate 

tr (e-iHoT/e+cPtT^HoT/e+cPoT (gl7) 

As before, one first has to check that the action of the semigroups can be substituted 
by the free evolution, 

\\ e -iHoT/e+cP,T p e e iH T/e+cP,T Q ^ _ Q -iH T / e / 1 Qp[r ] ^ 

< ^HoT/e^HoT/e+cP.T^HoT/e+cP.T^iHoT/e _ ^ . || O p[0 e ] || fl(w) 

< 2\c\e^ [ \\P^ H ^ +cP * s pl\\ tr dS • || V n\ B(H) . 

Jo 

The operator norm in the last line is bounded uniformly by the Calderon - Vaillancourt 
theorem. For the integral, Lemma [8.11 and the uniform boundedness of ||p lltr allow an 
application of dominated convergence, 

-T r-T 



limsup / ||P^e-^ oS/e+cP ^|| tr d5 = limsup / ||P^r iHo5 / e polltrd£ 

e^O Jo e^O Jo 

= limsup / \\p e e iHoS/€ (P\\dS = 0. 

e->0 Jo 



(8.19) 



The last equality is seen as follows: Since the spectrum of H is absolutely continuous, 
we first can consider a state with bounded spectral density, i.e. 

%P < C 2 (8.20) 
dE 

for some C > 0. Then for any ip e 7i, g(t) = ^,e lHot (p^ is the Fourier transform of an 
L 2 (IR) function, satisfying 

IMU2<v^c|M|. (8-21) 
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Therefore, the Cauchy-Schwarz inequality yields 

I | (ij,e iHoS/ ^)\dS < y/2iifC\\i/j\\y/€. (8.22) 
Jo 

Now we can write 

with {jpn) an e-dependent orthonormal basis and 2~^neN \ a n\ — K uniformly in e. Then 
(18T22|) implies ^ 

rT /-T 



/ ||Po^ o5A 0l|dS<^|<| / |<C,e^0)|dS 



(8.24) 



< kV^TC^ ^ (e -> 0). 



Approximating general <fi with states 0c of bounded spectral density, one can use the 
uniform boundedness of the operators p e lH ° s ^ e to show 

limsup / \\p e e iHoS/e <f>\\ dS < TK\\(f> - C ||, (8.25) 

e^O Jo 

where the right hand side can be made arbitrarily small. Thus one can in fact replace 
the semigroups by the unitary group, and then apply (18. 6p . 



limtr ( e -iH T/e + cP+T p e e iHoT/e+cP*T Q p[ ^ 

limtr (e- iHoT/t p e e iHoT/t Op[fl e ]) 



limtr (poOp[?7t]) 



= lim / W [p Q ] (x,p) 9(ex + pT,p)dxdp 

J^ 2d (8.26) 

= lim I e- d W [p £ ] (X/e, P) 6{X + PT, P)dXdP 

= [ g(X,P)6(X + PT,P)dXdP 

= / g(X - PT,P)9(X,P)dXdP. 
Adding (I8.16P and (I8.26p . one obtains the limit for general initial conditions, 
lim I f e ( y X,P,T)9( y X,P)dXdP 

= / g(X - PT,P)9(X,P)dXdP 



+ 2\r\ I [ 9(PS,P)\0(P)\ 2 dPdS 



8.27) 



V'.V 



A quantum particle source 
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The distribution in the last line is given by 

f°(X,P,T) =g(X-PT,P) + 2|c| I S(X - PS)\4>(P)\ 2 dS. (8.28) 

Jo 

It remains to show (13.191) . By (13.171) . the initial value is 

f°(X, P, 0) = lime'^ [pl\ (X/e, P) = g(X, P). (8.29) 

e— >0 

By testing with an arbitrary 9 E T> (li x x Mp) , one can check that f° also solves the 
differential equation, 

= ^ (J g(X, P)6{X + PT, P)dXdP^j 

+ ± (2 1 c\ £ J J (P(T - S) , P) 1 0(P) 1 2 dP dS^j 
= [ g(X,P)P-X7 x 0(X + PT,P)dXdP 



2\c\ [ [ P -V X 9(P(T - S),P)\4>(P)\ 2 dPdS 

JO JR d 

+ 2\c\ [ 9(0,P)\4>(P)\ 2 dP 



(8.30) 



(f(T),P-V x 9) v , v + (2\c\5(X)\<f ) (P)\ 2 ,9) i)i ^ 
-P-V x f°(T) + 2\c\6(XMP)\ 2 ,9) . 



□ 
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